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Abstract. The purpose of this article is to describe the two approaches to compute exact formulas 
(which are amenable to asymptotic analysis) for the probability distribution of the current of particles 
past a given site in the asymmetric simple exclusion process (ASEP) with step initial data. The 
first approach is via a variant of the coordinate Bethe ansatz and was developed in work of Tracy 
and Widom in 2008-2009, while the second approach is via a rigorous version of the replica trick 
and was developed in work of Borodin, Sasamoto and the author in 2012. 



1. Introduction 

Exact formulas in probabilistic systems are exceedingly important, and when a new one is dis- 
covered, it is worth paying attention. This is a lesson that I first learned in relation to the work of 
Tracy and Widom on the asymmetric simple exclusion process (ASEP) and through my subsequent 
work on the Kardar-Parisi-Zhang (KPZ) equation. New formulas can enable asymptotic analysis 
and uncover novel (and universal) limit laws. Comparing new formulas to those already known can 
help lead to the realization that certain structures or connections exist between disparate areas of 
study (or at least can suggest such a possibility and provide a guidepost). 

The purpose of this article is to describe the synthesis of exact formulas for ASEP. There are 
presently two approaches to compute the current distribution for ASEP on Z with step initial 
condition. The first (called here the coordinate approach) is due to Tracy and Widom j26l [271 ES] in 
a series of three papers from 2008-2009, while the second (called here the duality approach) is due 
to Borodin, Sasamoto and the author [5] in 2012. 

The duality approach is parallel to an approach (also developed in [5]) to study current distri- 
bution for another particle system, called g-TASEP. Via a limit transition, the duality approach 
becomes the replica trick for directed polymers. In fact, ASEP and g-TASEP should be consid- 
ered as integrable discrete regularizations of the directed polymer model in which the replica trick 
(famous for being non-rigorous) becomes mathematically rigorous. Underlying the solvability of 
g-TASEP and directed polymers is an integrable structure recently discovered by Borodin and the 
author [4] called Macdonald processes (which in turn is based on the integrable system surrounding 
Macdonald symmetric polynomials). It is not presently understood where ASEP could fit into this 
structure, but the fact that the duality approach applies in parallel for ASEP and g-TASEP compels 
one to look for a higher structure which encompasses both. 



2. Current distribution for ASEP 

ASEP is an interacting particle system introduced by Spitzer [2^ in 1970 (though arising earlier 
in biology in the work of MacDonald, Gibbs and Pipkin [18] in 1968). Since then it has become 
a central object of study in interacting particle systems and non-equilibrium statistical mechanics. 
Each site of the lattice Z may be inhabited by at most one particle. Each particle attempts to jump 
left at rate q and right at rate p {p + q = I), except that jumps which would violate the "one particle 
per site rule" are suppressed. We will assume q > p and for later use call q — p = "f and p/q = t 
(note that 7 > and r < 1). 
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There are two ways of constructing ASEP as a Markov process. The "occupation process" keeps 
track of whether each site in Z is occupied or unoccupied. The state space is y = {0, 1}^ and for 
a state rj = {r/^j^ez ^ Y , Vx = ^ there is a particle at x and otherwise. This Markov process is 
denoted r]{t). 

The "coordinate process" keeps track of the location of each particle. Assume there are only 
k particles in the system, then the state space = {xi < ■ ■ ■ < Xk} C and for a state 
X = {xi < . . . < Xk} G Xk the value of Xj is the location of particle j. We call Xk a Weyl chamber. 
Because particles cannot hop over each other, the ASEP dynamics preserve particle ordering. This 
Markov process is denoted x{t). 

In this article, we will be concerned with the "step" initial condition for ASEP in which every 
positive integer site is initially occupied and every other site is initially unoccupied. In terms of 
the occupation process this corresponds to having r]x{0) = lx>o (here and throughout 1_b is the 
indicator function for event E). Let Nx{r]) = '^y^^i]y and note that NQ{r]{t)) records the number 
of particles of ASEP which, at time t are to the left of, or at the origin - that is to say, it is the net 
current of particles to pass the bond and 1 in time t. 

Theorem 2.1. For ASEP with step initial condition and q > p, 



where FcuEis) is the GUE Tracy-Widom distribution. 

Remark 2.2. The distribution function -Fgue('S) can be defined via a Fredholm determinant as 

-?^gue(s) = det(/ - KAi)L2{s,oo) 
where Airy kernel K^i acts on L^(s, oo) with integral kernel 



For q = 1 and p = this result was proved in 1999 by Johansson [13] and for general q > p 
it was proved by Tracy and Widom [261 123 [SSj in 2009, and then reproved via a new formula by 
Borodin, Sasamoto and the author [5] in 2012. This result confirms that for all q > p ASEP is in 
the Kardar-Parisi-Zhang universality class |15j (see also the review |6]). 

In order to prove an asymptotic result (such as above), it is very useful to have a pre-asymptotic 
(finite t) formula to analyze. If the formula does not increase in complexity as t goes to infinity, there 
is hope to compute its asymptotics. Presently, there are two approaches to computing manageable 
formulas for the distribution of Noit). 



In [26] . Tracy and Widom start by considering the ASEP coordinate process x{t) with only k 

particles. In 1997, Schiitz [22] computed the transition probabilities (i.e.. Green's function) for ASEP 
with k = 2 particles. The first step in [26] is a generalization to arbitrary k. Let Py{x; t) represent 
the probability that in time t, a particle configuration y will transition to a second configuration x. 
As long as p 7^ 0, it was proved in |26j that 






3. The coordinate approach 
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where the contour of integration is a circle centered at zero with radius so small as to not contain 
any poles of Aa-. Here e(^) = p(,~^ + ~ 1 and 

Af^ = TT {SaB ■ {a, /3} is an inversion in a} , S^r = _ ^ ^'^"^/^ 

This result is proved by showing that that Py{x; t) solves the master equation for A;-particle ASEP 

^u{x; t) = L^u{x] t), u{x; 0) = 



Here L'^ is the generator of the /c-particle ASEP coordinate process. For k = 1, acts on function 
/ : Z ^ M as 

{L^f){x) = p [fix - 1) - fix)] + q [fix + 1) - fix)] . 

For k > 1, the generator L'^ depends on the location of vecx in the Weyl chamber, reflecting the 
fact that certain particle jumps are not allowed near the boundary of the Weyl chamber. 
Quoting a footnote in [26] : 

The idea in Bethe Ansatz (see, e.g. [HI [251 ED])) applied to one-dimensional k- 
particle quantum mechanical problems, is to represent the wave function as a linear 
combination of free particle eigenstates and to incorporate the effect of the potential 
as a set of — 1 boundary conditions. The remarkable feature of models amend- 
able to Bethe Ansatz is that the boundary conditions for A; > 3 introduce no more 
new conditions... The application of Bethe Ansatz to the evolution equation (mas- 
ter equation) describing ASEP begins with Gwa and Spohn [9] with subsequential 
developments by Schiitz [22] . 
To see this in practice, assume that one wants to solve 

— it(x; t) = L^uix] t), uix; 0) = uoix) 
for X in the Weyl chamber X^. 

Proposition 3.1. If v :'L^ x M solves the 'free evolution equation with boundary condition": 

(1) iA^;t) = Y:]=AL^'^A^-t),; 

(2) For all x G Z'^ such that x^+i = Xj + 1 for some 1 < J < A; — 1, 

pvixi,.. .,Xj,Xj+i -l,...,Xk;t) + qvixi,.. .,Xj + . . ,Xfc;t) - vix;t) = 0, 

(3) For all x E Xk, f (x; 0) = uoix), 

then for all t > and x G X^, ti(x;t) = vix;t). 

In (1) above, [L^]j means to apply the generator in the xj variable. In fact, some growth 
conditions must be imposed to ensure that u and v match (see Propositions 4.9 and 4.10 of [5]) but 
we will not dwell on this presently. 

This reformulation of the master equation involves only k — 1 boundary conditions and is amend- 
able to Bethe Ansatz - hence one is led to postulate equation ([1]). It remains to check the ansatz 
(i.e., that Pjjix; t) solves the reformulated equation). The A^j are just right to enforce the boundary 
condition. The only challenge (which requires an involved residue calculation) is to check the initial 
data, since there are a total of k\ integrals. 

The transition probabilities for /c-particle ASEP is only the first step towards Theorem 12. 11 The 
next step is to integrate out the locations of all but one particle, so as to compute the transition 
probability for a given particle Xm- The formula for the location of the m^^ particle at time t involves 
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a summation (indexed by certain subsets of {1, ... , A:}) of contour integrals. These formulas are a 
result of significant residue calculations and combinatorics. 

At this point we are only considering k particles, whereas for the asymptotic problem we want 
to consider step initial conditions. This is achieved by taking yj = j for 1 < j < k and taking k 
to infinity. After further manipulations, the m^^ particle location distribution formula has a clear 
limit as k goes to infinity. This is the first formula for step initial condition and it is given by an 
infinite series of contour integrals. 

In [27 { this infinite series is recognized as equal to a transform of a Fredholm determinant. By the 
simple relationship between the location of the m*'* particle of ASEP and N(j{t) (defined earlier), 
this shows that 

ITTI J [(,;T)m+l 

where the integral in ^ is over a contour enclosing (" = for < k < m — 1 and (a; r)„ = 
(1 — a)(l — to) • • • (1 — T"~^a). Here det(/ — C^i) is the Fredholm determinant with the kernel of 
K : L^Cr) ^ L^{Cr) given by 

p + qii' - i 

and the contour Cr a sufficiently large circle centered at zero. 

There remains, however, a significant challenge to proving Theorem 12. II from the above formula. 
As m increases, the kernel Ki has no clear limit, and the denominator term (("; T)m+i behaves widely 
as C, varies on its contour of integration. Much of |28j is devoted to reworking the above formula into 
one for which asymptotics can be performed. This is done through significant functional analysis. 
The final formula, from which Theorem 12. II is proved by asymptotics is that (leaving off the contours 
of integration) 

P(iVo(t) >m)= [ ^(/u; r)oo det(/ + fiJ), (3) 
where the kernel of J is given by 



k=—oo 



At,m,x(C) = -Xl0g(l-C) + Y3^+"ilogC• 



4. The duality approach 

Duality is a powerful tool in the study of Markov processes. It reveals hidden structures and 
symmetries of the process, as well as leads to non-trivial systems of ODEs which expectations of 
certain observables satisfy. In 1997, Schiitz [23] observed that ASEP is self-dual (in a sense which 
will be made clear below). In 2012, Borodin, Sasamoto and the author [5] used this duality, along 
with an ansatz for solving the duality ODEs (which was inspired by the work of Borodin and the 
author on Macdonald processes [4] ) to derive two different formulas for the probability distribution of 
A'^o (t) . The first was new and readily amendable to asymptotic analysis necessary to prove Theorem 
12. 1^ while the second was equivalent to Tracy and Widom's formula ([2]). 
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To define the general concept of duality, consider two Markov processes rj(t) with state space 
Y and x{t) with state space X (for the moment we think of these as arbitrary, though after the 
definition of duality, we will take these as before). Let E'' and E^' represent the expectation of these 
two processes (respectively) started from 7^(0) = rj and x(0) = x. Then r](t) and x{t) are dual with 
respect to a function H -.Y x X if for a\l r] £Y, x X and t > 0, 

[H{v{t),x)]=E-[H{7^,xm- 

One immediate consequence of duality is that if we define Ur,{x;t) to be the expectations written 
above, then 

—Ur,{x;t) = L*Ur,{x;t) 

where L* is the adjoint of the generator of x{t) and where the initial data is given by n^(x;0) = 
H{r],x). 

Schiitz [23] observed that if r]{t) is the ASEP occupation process and x{t) is the fc-particle ASEP 
coordinate process with p and q switched from the earlier definition, then these two Markov processes 
are dual with respect to 

i/(r?,f) = nT^-.-(^)ry,,. 

The adjoint of the generator of the p,q reversed x{t) process is simply L^. Schiitz demonstrated 
this duality in terms of a spin-chain encoding of ASEP by using a commutation relation along with 
the Uq[SU{2)] symmetry of the chain. A direct proof can also be given in terms of the language 
of Markov processes [5]. When p = q, t = 1 and this duality reduces to the classical duality of 
correlation functions for the symmetric simple exclusion process (see |17] Chapter 8, Theorem 1.1). 

As before, we focus on step initial condition, so that r]x = lx>i- Duality implies that Ugtepi^', t) ■= 
E^ [i7(7y(i),x)] solves 

d ^ 

— Ustcp{x;t) = L''Ustep{x;t), Ustep{x;0) = lxi>l JJt^"~^- (4) 

i=l 

The above system is solved by 

^fc(fc-l)/2 I- I- Zb ^ 

^stcp(x; t) = ... II — — — n hxjAzj)dzj, (5) 

where 

l + z 1 , z{p-q)^ 



KAz) = e^'(")* 



xAz) = e ^ ' — — — ■ — , e 



1 + z/t) t + z' {l + z){p + qz)' 

and where the contour of integration for each zj is a circle around — r, so small as to not contain 
— 1. In order to see this, we use the reformulation of the system in terms of the free evolution 
equation with boundary condition with ASEP given earlier in Proposition 13.11 Condition (1) is 
trivially checked since for each z, ^hx^t{z) = L^hx^t{z)- Condition (3) is checked via a simple 
residue calculation. Condition (2) reveals the purpose of the ^^Z^^b ^PPlyii^g the boundary 

condition to the integrand above brings out a factor of Zj — tzj. This cancels the corresponding 
term in the denominator and the resulting integral is simultaneous symmetry and antisymmetry in 
Zj and Zj^i. Hence the integral must equal zero, which is the desired boundary condition (2). 

The inspiration for this simple solution to the system of ODEs came from analogous formulas 
which solve free evolution equations with boundary condition for various versions of the delta Bose 
gas (see Section [5] for a brief discussion) . For the delta Bose gas and certain integrable discrete 
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regularizations, the formulas arose directly from the structure of Macdonald processes [4j. ASEP 
does not fit into that structure, but the existence of similar formulas suggests the possibility of a 
yet higher structure. 

A change of variables reveals some similarities to the integrand in ([1]) . Letting 



we have that 



ZA - ZB 



ZA - TZB P + qCaCb - 



1 + 

l + Zj/T 



hjcj,t{zj)dzj 



(6) 







The system @ could also be solved via Tracy and Widom's formula [T] for the Green's function 
for (see [12]) but the resulting formula would involve the sum of kl fc-fold contour integrals. 
Symmetrizing (jS]) via combinatorial identities, and making the above change of variables, one does 
recover that formula. The reversal of this procedure is a rather unnatural anti-symmetrization, 
which explains why ([5]) was not previously known. 

A suitable summation of H{r],x) over x gives t^^^M . Using this, and formula ©[Sj proves that 
for ASEP with step initial condition. 



E 



(27ri)fc 



n 



l<A<B<k 



ZA - ZB 
ZA - TZB 



n 



(7) 



where No^t) = NQ{r]{t)) and where the contour of integration for Zj includes 0, — r but not —1 or 
r times the contours for zj^i through Zk- This is to say, that the contours of integration respect a 
certain nesting structure. 

At this point the utility of having a single fc-fold nested contour integral formula for the moments 
of T^o(t) becomes clear. There are two ways to deform the contours of integration in ([7]) so as to all 
coincide with each other. The first involves expanding them all to be a circle containing — r and 0, 
but not —1. There are many poles encountered in the course of this deformation and the residues 
can be indexed by a partition. This leads to 



E 



kNo{t) 



Ahfc 
=1'"12" 



(1 



mi\m2l 





-1 


/.../det 


WiT^i — Wj 





where k-^l = (r; T)fc(l — r) is the r-deformed factorial, and A = (Ai > A2 > 



(8) 

> 0) is a partition 



of k (i.e. X] Aj = A;) with £{X) nonzero parts, and multiplicity rrij of the value j. The structure of 
these residues is very similar to the string states indexing the eigenfunctions of the attractive delta 
Bose gas (see Section [5]). 

The final step in the duality approach is to use these moment formulas to recover the distribution 
of Nolt). This is done via the r-deformed Laplace transform Hahn [10] introduced in 1949. The 
left-hand side of the below equation is the transform of r^"^*^ with spectral variable C- 



E 



00 

E 

fc=0 



^fe]£j^fcAfo(t)j 



(9) 



The right-hand side above comes from the left-hand side by expanding the r-deformed exponential 
inside the expectation (using the r-deformed Binomial theorem) and then interchanging the sum- 
mation over k with the expectation. This interchange of summation and integration is justified here 
for small enough because |r*^^°(*)| < 1 deterministically (in contrast to (I14p Section [5]). 
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Substituting ([8|) into the series on the right-hand side of ([9]) one recognizes a Fredholm determi- 
nant. The kernel of the determinant can be rewritten using a Melhn-Barnes integral representation 
and the result is (leaving off the contours of integration) 



E 



1 



(CT^0W;r), 



det(I + i^^), 



(10) 



where the kernel of is 

K(^{w, w') = 



1 

27ri 



TT 



-(-.) 



ds 



■yt-t- 



sin(— vrsj ~ ' g[T'^w) w' — T^w^ 

The r-Laplace transform can easily be inverted to give the distribution of No(t) and asymptotics 
of the above formula are readily performed (see Section 9 of [5j) resulting in Theorem 12.11 

There is a second choice for how to deform the nested contours in ([8]) to all coincide. The 
terminal contour of this deformation is a small circle around — r, and again there are certain poles 
encountered during the deformation. The combinatorics of the residues here is simpler than in the 
first case, and one finds the following Fredholm determinant formula 



E 



1 



where the kernel of K2 is 



.(Cr^oW;r), 



K2{w,w') 



det(/ - CK2) 

(C;7")oo 



(11) 



TW — W 

Performing the change of variables ^ and inverting this r-Laplace transform, one recovers Tracy 
and Widom's formula ([2]). As in Tracy and Widom's work, this formula is not yet suitable for 
asymptotics and must be manipulated significantly to get to the form of ([3]). 



5. Duality approach as a rigorous replica trick 

Besides the Schiitz duality, Borodin, Sasamoto and the author discovered that ASEP is also self 
dual with respect to 

k 

For k = 1 this shows that E[t^^(''(*))] solves the heat equation with generator L^. In fact, this is 
essentially Gartner's 1988 observation [8] that r^^^^^*)) solves a certain discrete multiplicative sto- 
chastic heat equation. A multiplicative stochastic heat equation has a Feynman-Kac representation 
which shows that the solution can be interpreted as a partition function for a directed polymer in 
a disorder given by the noise of the stochastic heat equation. 

In 1997 Bertini and Giacomin [2j showed that under a certain "weakly asymmetric" scaling 
^N^{r){t)) converges to the solution to the continuum multiplicative stochastic heat equation (SHE) 
with space-time white noise ^(x,t): 

j^Z{x, t) = l^Z{x, t) + Z{x, t)e(x, t). 

This convergence result did not include when r/(0) is step initial condition and was extended to 
that case by Amir, Quastel and the author [ij. The corresponding initial data for the SHE is 
Z{x, 0) = 6x=o where S is the Dirac delta function. 
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Duality of ASEP translates into the fact that the moments of the SHE solve the attractive one- 
dimensional imaginary-time delta Bose gas (Lieb-Liniger model with delta interaction) [14J . Define 
Z(x; t) = E [Z{xi,t) ■ ■ ■ Z{xk,t)] for Z with Sx=o initial data. Then Z solves the system 

d - - - ^ 

-Z{x; t) = HiZ{x; t), Z{x; 0) = J] '^-.=0, (12) 

i=i 

where Hf^ is the Lieb-Liniger Hamiltonian with delta interaction with strength k € M: 

j=l 3 Kj 

The Lieb-Liniger model with delta interaction was the second system solved by the Bethe Ansatz 
(over 30 years after Bethe [3j solved the spin- 1/2 isotropic Heisenberg model). This was accomplished 
by Lieb and Liniger in 1963 for the repulsive system (in which Hi is replaced by H,^ for k < 0). 
A year later, McGuire similarly solved the attractive system (k > 0). In their context, solving the 
system meant writing down eigenfunctions for ff^- The structure of the eigenfunctions for k < 
and K > are different. In the attractive k > case there are extra eigenfunctions which are called 
"string states" due to the strings of quasi-momenta with which they are indexed. Completeness of 
these eigenfunctions was not shown until later [TJ [TTl [29l [191 120] ■ 

For the purposes of understanding the moments of the SHE it is not necessary to diagonalize Hf^, 
but rather just to solve the system (fT2]) for n = 1. Just as with ASEP, this system can be written 
as a "free evolution equation with boundary condition". The free evolution is just according to the 
k variable Laplacian and the boundary condition is that for all 1 < j < /c — 1, 

d d \ ,^ .1 
; K v(x; t)\ ^ =0. 

JXj dXj+i J ^ '\xj^Xj+i 

This system can be solved via an analogous formula to For xi < . . . < xi^. and k G M, 

^ ^ l<A<B<k ^ ^ j=l I ) 

where the contour of integration for Zj is along aj + iR for any ai > 02 + k > ■ ■ ■ + {k — 1)k. 
When K < all the aj can be chosen as and hence the integral occurs on iM, whereas for k > 
the contours must be spaced horizontally. In the k > case, the contours can be deformed to iM. 
The singularities and associated residues encountered have a very similar structure to those seen 
earlier in ([8]) in the context of the first ASEP contour deformation. The disparity between residue 
combinatorics accounts for the difference in the structure of the eigenfunctions and the occurrence 
of string stated for At > 0. In fact, Heckeman and Opdam's 1997 proof of the completeness of the 
Bethe Ansatz relied on a formula equivalent to ()13p . 

Given expressions for all of the moments of the SHE, one wants to recover the distribution of 
Z{x,t). Since Z(x,t) is nonnegative, its Laplace transform characterizes its distribution. Naively 
one writes 



E 



^ C'E[Z(t,x)fe] ^ ^^^^ 
fc=0 ^' 

However, the right-hand side is known to make no mathematical sense and the interchange of 
expectation and summation is totally unjustifiable. The moments of the SHE grow like e'^^^ and 
thus the right-hand side is extremely divergent. One can see that cutting off the summation also 
fails to remedy the situation in any way. 
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What should be clear now is that ASEP is an integrable discrete regularization of the SHE (or 
equivalently the KPZ equation) and the duality approach to solving it is a rigorous version of 
the replica trick for the SHE. By taking the weakly asymmetric limit of the r-deformed Laplace 
transform formulas described above, one finds a Fredholm determinant formula for E |^e''^(*'^)] . This 
can be done from either the new formula (llOp in [5] or Tracy and Widom's formula It appears 
that (llOp is very amendable to asymptotic analysis. 

Using ([3]), the derivation of the Laplace transform of Z{t,x) involves extremely careful asymp- 
totic analysis which was performed in 2010 rigorously by Amir, Quastel and the author [IJ and 
independently and in parallel (though non-rigorously) by Sasamoto and Spohn [21j. Very soon 
afterwards, Calabrese, Le Doussal and Rosso, as well as Dotsenko showed how to formally recover 
this Fredholm determinant from summing the divergent series on the right-hand side of (jl4p . The 
formal manipulations of divergent series that goes into this can be see as shadows of the rigorous 
duality approach explained above for ASEP. It can also be seen as a shadow of a parallel duality 
approach for q-TASEP [H[5], another integrable discrete regularization of the SHE. 
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